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Abstract 
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fer matrix of the sl(3) chiral Potts model for q 2 = —1. For the homoge- 
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dition, a direct nested Bethe ansatz has also been developed for this 
case. 
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1. Introduction 



The discovery of the chiral Potts model (CPM) and new iV-state solutions 
of the Yang-Baxter relation [1-4] is one of the most impressive results in 
the theory of two-dimensional lattice integrable systems [5]. A remarkable 
feature of this model is that the spectral parameters belong to a high genus 
algebraic curve. The reason of this fact remained unclear until it was shown 
that the weights of the chiral Potts model naturally appear as intertwiners of 
the cyclic representations for si (2) L-operators [6] related to the six- vertex 
model [7]. 

The next step was to generalise the chiral Potts model for the sl(n) case 
for q N = —1 [8]. As a result, a new A rn_1 -state family of generalised chiral 
Potts models has been obtained with Boltzmann weights satisfying the Yang- 
Baxter relation [9, 10]. 

Then Bazhanov and Baxter [11] made a remarkable observation that the 
sl(n) chiral Potts model with N = 2 is related to the integrable three- 
dimensional Zamolodchikov model [12,13]. Therefore, a new link between 
two quite different regions of the theory of integrable models has been re- 
vealed. It is well known that the tetrahedron equations provide sufficient 
conditions for a commutativity of the transfer matrices in three dimensions 
[14,15]. As a result, a new iV-state three-dimensional family of integrable 
models has been discovered with Boltzmann weights satisfying the tetrahe- 
dron equation [16, 17]. 

The Zamolodchikov model attracted a lot of attention even before the 
appearance of the chiral Potts model. Its partition function in the thermo- 
dynamic limit has been calculated in [18] and the related hamiltonian has 
been studied in [19]. It was shown in [19] that the Zamolodchikov model with 
two layers is equivalent to the critical two-dimensional free-fermion model. 
However, the structure of the spectrum of Zamolodchikov model with three 
and more layers remains unclear. 

In this paper we try to take several steps towards the understanding of 
the structure of the spectrum of si (3) chiral Potts model which is equivalent 
to the "modified" (in the terminology of [11, 18]) three-layer Zamolodchikov 
model. Starting from the general case we carefully analyse the homogeneous 
case of the sl(3) chiral Potts model where a usual nested Bethe ansatz tech- 
nique [20] can be applied. 

The paper is organised as follows. In Section 2 we give a full description 
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of the s/(3) chiral Potts model and introduce all necessary notations. In 
Section 3 we describe a "fusion" procedure [21-23] for the s/(3) chiral Potts 
model for q 2 = — 1 (N — 2) and obtain the closed functional equations for 
eigenvalues of the transfer matrix. Section 4 is devoted to the spectrum of 
the model in the completely homogeneous case. In Section 5 we apply the 
nested Bethe ansatz technique to obtain eigenvalues of the transfer matrix 
for the homogeneous case. In the last section we discuss our results and 
further directions for investigation. 



2. The 5/(3) chiral Potts model 

In this section we give all necessary definitions for the s/(3) chiral Potts 
model with fixed N. All definitions and formulas of this section are just a 
specification of the sl(n) case for n = 3 [10]. 

Consider an oriented square lattice C and its medial lattice CJ (shown 
in Figure 1 by solid and dashed lines, respectively). The oriented vertical 
(horizontal) lines of CJ carry rapidity variables Pi,p% {t\, t 2 ) in an alternating 
order (note that the orientations of rapidity lines shown by open arrows 
alternate too) . The edges of the lattice C are oriented in such a way that all 
the NW-SE edges have the same (NW-SE) direction while the NE-SW edges 
are oriented in a checkerboard order. 

Each rapidity variable is represented by three 2- vectors (hf(p), h~(p)), 
i = 1,2,3 which specify the point p of the algebraic curve T defined by 
relations 

(ht(p) N )- K ..(K(P) N \ v, 7-123 (21) 
\K(P) N ) ~ t3 \hj(p) N )' V ^- i ' 2 > d - V- L ) 

where Kij are 2x2 complex matrices of moduli satisfying relations 

detKij = 1, K u = K tJ K jk K kt = 1. (2.2) 

Hereafter we imply that indices k, . . . run the values 1, 2, 3 modulo 3. 
Below we need the automorphism r on the curve T defined as follows 

h+(r(p)) = h+(p), hj(T(p))=uhj(p), j = 1,2,3. (2.3) 

On each site of the lattice C place x Zn spins, which are described 
by a local variable 

a = (ai, a 2 ), ol\, a 2 — 0, 1, . . . , N — 1. (2.4) 



3 



There are two kinds of neighbouring local state pairs depending on the 
relative orientation of the dashed and solid lines as indicated in Figure 2, 
with states a and (3, and Boltzmann weights W pq (a,(3) and (W qp (a, 
on the edges of C. 

The function W pq (a,(3), a, (3 E Z N x Z N is defined by the following 
relations 

W m (a,(3) = uj Q ^g pq (0,a-(3), u = exp(— ), (2.5) 

where 

Q(a, (3) = - ai) + (3 2 {f3 l - on + (3 2 - a 2 ), a, (3 E Z N x Z N (2.6) 
and the function g pq (0, a) has the following form 

ai+»2— 1 

II (h+(p)h^(q)-h+(q)h^(p)u-P) 
9 P M a) = (2.7) 

i=l ft=0 

We choose a normalisation of W pq (a,(3) as 

W M (0,0) = 1. (2.8) 

Then it is easy to see that 

W ppM) ^, ^ S {J; 1=^^^ (2.9) 

The function W P9 (ai, /3) satisfies the inversion relation 

£ W pg (a,/3)TF gp (/3, 7 ) = 5 Q , 7 $ P9 , (2.10) 

where 

Xp~ t=\ X i\P) ~ Xi{q) N K(p) £=1 

(2.11) 
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Now let us suppose that our lattice C has M sites in a horizontal direction 
and L sites in a vertical one (L should be even). As usual we imply cyclic 
boundary conditions in both directions. 

Let us denote the spin variables of three consecutive rows as <7i, . . . , <Jm, 
a[,...,a' M and a", . . . , s" M (see Figure 1). Then we can define two N 2M xN 2M 
transfer matrices T pi and T P2 of the length M 



M 



Pklfcfe = U W PirA<7i,<7'i)Wri Pl (<7>i + l), (2-12) 



i=l 



\r i<-^ _ fx Wp2rJ^i±v»^ , 9 , o\ 

1=1 " P2V1 \ u l J u iJ 

Let us fix the rapidity variables r±, r 2 and use more simple notations T pi and 



T 

P2- 



The partition function is 

Z = Tr(T Pl T P2 ) L / 2 . (2.14) 

In the next sections we derive several functional relations between T pi 
and T p2 . 

We can construct two /^-matrices from the weights W pq (a, (3). Define (see 
Figure 3) 

?fe[ffiVl|= «^«^ (2.15) 

This _R-matrix satisfies the Yang-Baxter equation 

Sn(Pi,P2; qi, q2)S 13 (pi,P2; r 1 ,r 2 )S 2 3(qi, <?2; ri, r 2 ) = 

= 5 , 23(gi,g2;^i,^)5'i3(pi,P2;n,^2)5'i2(pi,P2;gi,g2)- (2.16) 

It leads us to the following commutation relation for the product of the 
transfer matrices T p2 T pi 

Tp2T pi Tq 2 T qi — T q2 T qi T P2 T pi , (2.17) 

where p±, p 2 , qi, q2 are four points on the curve T. 
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Similarly one can introduce the second i?-matrix S (see [10], for example) 
and prove another commutation relation 

-^pi-^P2^~<?i-^<32 = ^9i^<22^pi^p2- (2-18) 

However, all commutation relations between the transfer matrices are the 
consequence of the only equation which is called "star-star" relation [24], [11]. 
To write it down define two "stars" (see Figures 4,5) 

W»» (a, (3, 7 , S) = £ ^(^^ 2 (7^)W r2Pl (a,^) (2 JQ) 

and 

w^MM,i) = J f^W^f«W™(s.°) , (2 . 20) 

In these notations the star-star relation can be written as (see Figure 6) 

Wp 2Pl (5, a)W r2ri (5, 7 )W^ (a, (3, 7 , 5) = 
W™ (a, (3, 7 , 5)W nri (a, (3)W P2Pl (7, (3). (2.21) 

The Yang-Baxter equation (2.16) can be easily obtained by repeated appli- 
cations of (2.21). 

Now using (2.21) we obtain one more "commutation" relation between 
the transfer matrices T Pl , T P2 . To do that define the shift operator 

M 

Kt:h - U^, <+1 (2.22) 

i=l 

Then using cyclic boundary conditions in a horizontal direction and relation 
(2.21) we obtain 

T pl T P2 U {d) X- 1 = U {d) X- l T P2 T Pl , (2.23) 
where is the diagonal matrix independent of pi, P2 

M 

P {d) \a\::::^ = W 5 ^w r2 ^a l+l ). (2.24) 
i=i 
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3. Functional relations at q 2 = — 1 



Here we consider only the case q 2 = uo = — 1. The product of the transfer 
matrices T Pl T P2 can be easily rewritten in terms of "star" weights (2.19) 

M 

T P1 T P2 = n W^(a u a i+1 , a> +1 , a[). (3.1) 
i=i 

We can specify the horizontal rapidities in such a way that the left and right 
"gauge" factors W Pl:P2 (a, (3) in (2.21) become degenerate. Namely, there are 
two simplest choices to do that 

P2 = r\p), p 1 =p, A = 0,1, (3.2) 

where r is defined in (2.3) and the case A = corresponds just to p 2 = P\ = p. 

First let us set p 2 = Pi = P- Then using (2.21) and the explicit form of 
W plP2 (a, (3) it is easy to see that we have the only non-zero matrix elements 
for W^ r2 (a, (3, 7, a) if (3 — 7. Similarly setting p 2 = t(pi) we obtain that the 

only non-zero matrix elements for W^f 2 {a, (3, (3, 5) are when a — 5. 

Therefore, we obtain for T p T T x^ that if 07 = a" for some I — 1, ... ,M 
then 

[T p T t xm]%;:::;% = > if for some J^I, J=1,...,M 

(3.3) 

As a result we can split T P T T \^ into "diagonal" and "non-diagonal" parts. 
To do this explicitly consider two 4x4 projectors 

Pfjs = \, P;,s = Ks-\i 7,5 = 1,. ..,4. (3.4) 

Now let us consider i?-matrix (2.15) as 4x4 matrix with respect to horizontal 
indices with fixed vertical ones and denote it as S piP2 (a, (3). Then we have 

P + S pp (a, f3)P + = P + S pp (a, (3), P~S pp (a, (3)P~ = S pp (a, (3)P~, (3.5) 

P + S pT(p) (a,f3)P + = S pT(p) (a,f3)P + , P'S^a, (3)P~ = P~ S pT(p) (a, (3). 

(3.6) 
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As a consequence we have the following decomposition of the product T p T T x^ 

T p T r , (p) =Tr(P + V, (p) (a,/3)P + ) M + Tr(p-V, (p) (a,/3)p-) M . (3.7) 

It is not difficult to check that 

Tr(P + S pT , (p) (a,(3)P + r = $^7, (3.8) 

where I being the identity matrix 4 M x 4 M , Q pq is defined in (2.11) and all 
subscripts should be considered modulo 2. Denote the second term in (3.7) 

as 

T w (p;r 1+X ,r 2+X ) = Tr{p-^ pT , {p) {a^)P~) M . (3.9) 

The transfer matrix (3.9) has nonzero matrix elements only if 07 7^ a", 
W = 1,...,M. 

It is not difficult to construct local "fused" weights for the transfer matrix 
TW(p- n + x, r 2+x ). First let us write the following decomposition for P~ 

^ 7 "> = !>(*> tMM), c(i,a) = lu iai+ ^ i - 1 ^, i = 1, 2, 3, a e Z 2 x Z 2 . 
i=i ^ 

(3.10) 

Now consider the family of "fused" L-operators which act in the tensor prod- 
uct of the auxiliary space C 3 and the quantum space C 4 

L ( ij\p;r 1 ,r 2 ) = A p - ru r 2 E c(i,7)c(j\ <S)^(p,T A (p);r 2 ,ri), (3.11) 

■y,SeZ 2 xZ 2 

A ™ = ' (r2) * (p) + 1 {p)t{r2)) 

(3.12) 

and 5 ' (pi,p 2 ; ^2, ri) is defined in (2.15). 



Introduce the transfer matrix t^ x \p; r\, r 2 ) constructed from the "fused" 
perators L^\p\ r 1? r 2 ) over the auxiliary space C 3 

M 

t (x \p; n, r 2 rk::A = E EII n> ( 3 - 13 ) 

i\ Jjk a=l 
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It is easy to check that two transfer matrices T^ x \p; r\, r 2 ) in (3.9) and 
t^ x \p;ri,r 2 ) are related as follows 

-TT7 — t (A) (p;n,r 2 ). (3.14) 

Therefore, the problem of calculating eigenvalues for the product of two 
transfer matrices T p T T x^ in (3.7) is reduced to a calculation of the eigenval- 
ues for the transfer matrix t^{p; r l5 r 2 ) in (3.13). 

Let us give explicit formulas for the matrix elements of L^\p; r 1? r 2 ). We 
have 

L^\p;n,r 2 ) =ut i . 1 (p;r 1 ,r 2 )Xi +1 + Ui_ 1 (p;r 1 ,r 2 )X i - 1 , (3.15) 

L i%i(p; r i, r 2 ) = Z iii+1 (^ ) 1 _ A (p; n, r 2 )X m _ A + «fc A) (p; n, r 2 )X i _ 1 ), 

(3-16) 

4?-i(p;n,r 2 ) = Zi ii _i(^i 1 A) A (p;ri,r 2 )Xi_i +A + w^ ) 1 (p;ri,r 2 )X i+ i), (3.17) 
where 

«f (p;ri,r 2 ) = ^ A *£(p)^(n)> (3.18) 

uPfan,^) = h?(p)ht(r 2 )h- +1 (n)ht +1 (p)h-_ i(p)hti(n) (3.19) 

4 ?) (p;ri,r 2 ) = ^(^^(r,)^^)^^)^!^)^.!^), (3-20) 

the 4x4 matrices X± and act in Z 2 x Z 2 and have the following matrix 
elements 

< a\X 1 1/3 > = 5 q1 ,/3i5q 2 ,/32-1) < a |^2|/3 > = £ai,/3i-l<5a2,/3 2 ; (3-21) 

< >= cc? a2 5 aiA <5 a2i/ 3 2 , < a|Z 2 |/5 >= cu ai (5 aiA 5 a2i/32 , (3.22) 

X 1 X 2 X 3 = 1, Z 3 = l, Z lj = Z t Zi 1 , i^j = 1,2,3. (3.23) 
They satisfy usual relations 

X,Z jfc = c^-^-fcXi, %^ fc Z fei = 1. (3.24) 

In fact, it is convenient to introduce the "gauge-transformed" L-operators 

L§\x p ; ri ,r 2 ) = [n^Wj^KiWl^^^rx.^feCp)]^, (3.25) 

i=l 
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where 

a?) = (Up), Up), Up)) = Mp), W(p))- (3.26) 

and x p , Xi{p) are defined in (2.12). 

It is easy to see that L-operators (3.25) depend on the rapidity variable 
p only via the combination x p . Therefore, we have 

t^(p; ri ,r 2 ) = (l[hi(p)) M t^(x p ,r u r 2 ), (3.27) 
i=i 

where t^(x p ; n, r 2 ) is the transfer matrix constructed from the L-operators 
(3.25) and obviously just a polynomial in the variable x p of the degree M. 

It follows that we can consider the parameters x p , hf(ri), hf(r 2 ) as in- 
dependent ones and do not take into account the surface equations (2.1). 

L-operators (3.25) satisfy the following Yang-Baxter equation 

E ^n,n(*/y)L^)^ty 2 )= E ^^)^^)B^^ ]3 {x/yl 

12, 32 «2J2 

(3.28) 

where we omit a dependence on r\ , r 2 in (x;r 1 ,r 2 ) and R il irjl j2 {x) coin- 
cides with the deformed trigonometric s/(3) R- matrix. Namely, let us assume 
for a moment that q is arbitrary and introduce i?-matrix [25] 

Ril,i2;ji,j2 i X i OS P) = ^iii2^jij2^iijii ( l i){x -\- X q ) + 
+ ^hi2^jij2Phji( X ~ X ) + ^Hj2^2jl(l ~~ ^iiji) a hi2{ X ), (3.29) 

where 

[0, i=j; 
Pa = Pijpji = 1, on = < (q - q^)x, i < j; (3.30) 

{ (q-q-^x' 1 , i>j. 

Then 

R{ &2-,h,ni x ) = RhMjuhi^^Q, P) (3-31) 

with 

q = i, p hJ+1 =p^ = l ^) x+ \ j = 1,2,3. (3.32) 

Now using a fusion technique for the i?-matrix (3.29) we can construct func- 
tional equations for the transfer matrix t^(x p ; r±, r 2 ) . Actually this proce- 
dure is quite similar to the fusion for the sl(3) trigonometric i?-matrix, so 
we do not give a detailed derivation of the functional equations. 
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We can define L-operators related to the "antisymmetric" representation 

~ (-l) A+5l ' 2 ^((-l) A ^ n, r 2 )LS((-l) 1+A x; n,r 2 )} , (3.33) 

where indices {i,k,l} and {j,m,n} are even permutations of {1,2,3} and 
we define 

0f (x; r) = ([] K(r) ±x]l h+(r)). (3.34) 
i=i i=i 

Let ri, r 2 ) be the transfer matrix constructed from the L-operators 

(3.33). 

Omitting a dependence of t^ x \x p ;ri,r 2 ), t^ X \x p ;ri,r 2 ) on r±, r 2 one can 
show that the following system of functional equations holds 

t^(x p )t {x \ux p ) = <f> (xp,r 1+x , 2+ x) + 

+fc(x p ,r 1+x ) M t (X \x p ) + <f>t(x p ;r 1+x ) M t {X \ux p ), (3.35) 
t W (x p )t W (ux p ) = <$>o(x p ;r 2+XA+x ) + 

+^(x p ;r 2+x ) M t^(x p ) +<Pt(x P ;r 2+x ) M t( x Xcux p ), (3.36) 

with 4>i{x p \ r) defined in (3.34) and 

$o(w',r) = Af + Af + A 3 M , (3.37) 
where Aj, i = 1,2,3 are three roots of the following cubic equation 

A 3 + a\ 2 + b\ + c = 0, (3.38) 

3 3 h + ( r '\2 3 3 h -( r l\2 

3 ht{rf "h-{rf 2 ^h-{T'f^h-{rf 3 /*+(r) 2 



=1 [^M]- 2 ^ (r) 2 P j=i ^') 2 W) 2 P ^ ^') 2j 

(3.40) 



3 



c=«n - n mo 2 )**; n w 2 - n mo 2 )) 2 - (3.41) 



1=1 1=1 1=1 1=1 
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A system of functional equations similar (3.35-3.36) for "factorized" si (3) 
L-operators has been obtained in [24] and discussed in [11]. However, for 
that case there is no Z 3 symmetry with respect to the indices of L-operators 
and as a result, the functional relations will have a slightly more complicated 
form and involve explicitly shift operators in the "quantum" space. 

It is easy to see that the system of functional equations (3.35-3.36) is 
invariant under the replacement T\ — > r 2 , t^ x \x p ;ri,r 2 ) — > t^\x p ; ri, r 2 ). 
However, it is not true in general that all eigenvalues of the transfer matri- 
ces t^(x p ), t^ X \x p ) satisfy t <yX \x p ;r 1 ,r 2 ) = t^(x p ;r 2 ,ri). The spectrum of 
t^(x p ), t^\x p ) involves also "non-symmetric" solutions of (3.35-3.36). 



4. The homogeneous case r\ = r<i. 

A great simplification occurs if we restrict ourselves to the case of the com- 
pletely homogeneous model r\ = r 2 = r. In this case all roots in (3.38) 
coincide and the function & (x p ;r,r) has a very simple form. Let us define 

-Si (41) 

It is easy to see that for this case the transfer matrices t(x p ; r, r) and t(x p ] r, r) 
are effectively the functions of a single variable x = x p /x r . Namely, 

3 3 

t(x p ;r,r) = (x r ]J h+(r)) M t(x p /x r ), t(x p ;r,r) = (x r J[ hj(r)) M t(x p /x r ). 
i=i i=i 

(4.2) 

Then omitting the superscript (A) the system of functional equations (3.35- 
3.36) can be rewritten as follows 

t(x)t{-x) = 3(1 - x 2 ) M + (1 - x) M t(x) + (1 + x) M t(-x), 
t(x)t(-x) = 3(1 - x 2 ) M + (1 - x) M t(x) + (1 + x) M t(-x), (4.3) 

where t(x) and t(x) are polynomials in x of the degree M. 
Let us define 

u(x) = t{x) + {1 + x) M , u(x) = t(x) + (1 + x) M . (4.4) 
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Then we have from (4.3) 

u(x)u(—x) = u(x)u(—x), (4.5) 

u(x)u(-x) = (1 - x) M (u(x) + u(x)) + (1 + x) M (u(-x) + u(-x)). (4.6) 
The solution of (4.5) can be written as follows 

k 

u(x) — a(x)b(x), u(x) — a(—x)b(x), a(x) — J\(v,i — x), k — 0,...,M, 

i=i 

(4.7) 

Ui, i — 1, . . . , k - the roots of the polynomial a(x) of the degree k and b(x) is 
the polynomial in x of the degree M — k. 
Then we have from (4.5-4.6) 

c(x)c(-x) = (1 - x 2 ) M (a(x) + a(-x)) 2 , (4.8) 

c(x) = a(x)a(-x)b(x) - (1 + x) M (a(x) + a(-x)), (4.9) 

where c(x) is the polynomial in x of the degree M + k. 

In principle it is straightforward to solve (4.8-4.9) in terms of Ui. First 
we need to find the roots of the polynomial a(x) + a(—x): 

a(x) + a(-x) = A J] (v 2 - x 2 ), (4.10) 
i=i 

where \_k/2\ is the integer part of k/2. 

Now a general solution of (4.8) has the following form 

h W2\ 

c{x) = \{l-x) n {l + x) M - n \{{v l -x) n (vi + x)x 

1=1 i=h+l 

h Lfc/2J 

x[](^-^) II fa + x), n = 0,...,M, < Zi < h < Lfc/2J. (4.11) 

i=l i=Z 2 +l 

The last step is to substitute (4.7, 4.10-4.11) into (4.9) and demand that b(x) 
should be a polynomial in x. It gives us a closed system of equations on Ui 
and V{. 
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As a result we come to the system of equations on Ui and vf. 



k «,--«.- k Lfc/2J 



II - = -1, i = 1, • • • , L^/2J , = II (^--«?), i = 1, • • • , A;, 

(4.12) 

[[— 11 = -(i ) , n = 0, ...,M,hJ 2 = 0,...,[k 2\. 

j=1 Vj +Ui j=h+1 Vj -Ui ^l-Ui' 

(4.13) 

It is not difficult to solve (4.12-4.13) for k — 0, 1 and to obtain some 
eigenvalues of the transfer matrix T p T p at r\ = r 2 for any M. However, in 
general the system (4.12-4.13) is a transcendental system of equations and 
can not be solved explicitly for any values of k. 

Therefore, (4.7-4.13) together with (4.1-4.2, 4.4) give the solution for 
the spectrum of the transfer matrix T p T p in the homogeneous case r\ = 
r 2 . However, we should say that this case is quite restrictive and, in fact, 
all dependence on the vertical rapidity r and matrices in (2.1) can be 
absorbed by a redefinition of the parameter x p . 



5. The Bet he ansatz technique for "fused" L- 
operators. 

In the previous section we considered the solution to the functional relations 
(3.35-3.36) for the case pi = p 2 = p, r\ = r 2 = r. However, a direct 
Bethe ansatz method can be also developed for this case. Let us make some 
transformations of the L-operators given by (3.15-3.17) (here we consider 
only the case A = 0). Define the //-operators: 

Llj = k Up)/Ur) L^(p;r,r) fc(r)/fc(p), (5.1) 
where is defined in (3.26), 

« = (h + (p)h-(p)h + (r)h-(r))-^, h^p) = f[ Kip). (5.2) 

i=i 

The transfer matrices for L^- and //-operators differ only by the scalar factor 
k m . Then all matrix elements of V are the functions of a single parameter 
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x = (xp/xr) 1 ^ 2 . To obtain the simplest form of L' let us make the following 
equivalence transformation: 



/I 



— C L'^ C 1 , 



C = 



1\ 



1 1 
1 1-1-1 

1-1 1-1 

\i -l -l 1/ 



(5.3) 



Then we have for the matrix elements of the L-operators 

Lii(a,(3) = {x}(5 afi 5p fi - SajSfri) + [x](5 a> i +1 5p ji+1 - (5.4) 
Lij(a, 13) = 2x eij (5 a , k 8p, - 5^5^), (5.5) 

where {x} = x+l/x, [x] = x—l/x, the indices i — in (5.4) are defined 

by the cyclic permutation of 1, 2, 3 while k in (5.5) is any permutation of 
1, 2, 3 and 

1, i < j\ 



-1, i > j. 



(5.6) 



In (5.4-5.5) we imply that the enumeration of rows and columns is given by 
(0,1,2,3). 

From (5.4-5.5) we can easily observe that i) the representation matrices 
of the diagonal L-operators are diagonal and for all of them the component 
(0,0) is equal to {x}; ii) the non-diagonal L-operators have zeros in the 
0-rows. 

Therefore, we can conclude that for the transfer matrix 



rp/3l,...,PM 

ai,...,(3 M 



]T L ili2 (a 1 ,p 1 )... («M,/3 



Mi 



(5.7) 



n,...,t M 



all matrix elements T^"'^ M = 0, for which there is at least one 0-component 
among the set a±, . . . , ocm and all indices jk run the values (1, 2, 3). So, the 
transfer matrix (5.7) has the block-down-triangular form. 
For example, the simplest block has a dimension 1: 



i 0,...,0 



3{x} M . 

Further there are M blocks 3 x 3 of the form: 

{x} M-l T W^ 



rpO,...,0,j k fl,...,0 

J 0,...,0,i fc ,0,...,0 



(5.8) 



(5.9) 
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where is the one-site transfer matrix for the "reduced" L-operators which 
can be obtained from L by removing all O-components, i.e. 

Lij(ii,ji) =L ij (i 1 ,j 1 ) = (1,2,3). (5.10) 

The next M(M — l)/2 blocks in the transfer matrix have the dimension 

3 2 x 3 2 : 

rp0,...,0,j kl ,0,...,j k2 ,0,...,0 _ f„\M-2 T (2)3kiJk2 (r, 

1 o,...,o,i kl ,o,...,i k2 ,o,...,o - \ x l 1 i kl ,i k2 - y°- Ll ) 

where T^> is the two-site transfer matrix for the "reduced" L-operators. In 
the n-th step we have blocks obtained from the n-site transfer matrices. 
In the last step we have one block: T^ M ^'[''^ . As a result we obtain a 
decomposition of the initial 4 M -dimensional space into the direct sum of 
the subspaces Y,n=o C'm3 t l . Therefore, the initial problem of calculation of 
the spectrum for the transfer matrix (5.7) has been reduced to the spectral 
problem for the transfer matrices T^ n \ 

It is easy to see that the L-operators given by (5.10) coincide with the 
i?-matrix (3.29) for the si (3) model. Hence, one can use the standard nested 
Bethe ansatz technique [20] for this model. We give here only the final 
result for the eigenvalues for the deformed sl(3) model with the arbitrary 
deformation parameters q and pij. 

am = j^L n «& + rn^ri <)n < > + %i\ «.<f-)i, 

Pl,3P2,3 i=l X Pl,2 i=i Vi k =l X Pl,2 k=l Zk 

(5.12) 

where < b < a < n, w(x) = W and the two sets of parameters yi, . . . , y a 
and Z\, . . . , Zf, should be defined from the system of the Bethe ansatz equa- 
tions: 

where p = Pi,2P2,3P3,i- However, we need only a particular case of this solu- 
tion 

Q = h Pjj+i = -Pj+i,j = ~h 3 = 1, 2, 3. (5.14) 
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In this case the formulas (5.13) become much more simple: 



w'^r = (-1)— 1 n «/(-), n ^'(-) = (-1)-- 6 - 1 , w(x) = m. 

fc=l ^* i=l f» L X J 

(5-15) 

In comparison with the equations (4.12-4.13) the equations (5.15) do not con- 
tain the redundant solutions. The reason for this can be easily understood. 
Namely, if we consider another type of //-operators which can be obtained 
from the initial ones by the transposition of the representation matrices, one 
can conclude that the transfer matrix for them satisfies the same functional 
relations (4.3). The resulting equations of the Bethe ansatz for this case can 
be obtained from (5.13) by the following substitution: 

x -> x~\ yi -> y' 1 , z k -> z^ 1 , p u -> p jti . (5.16) 

After this we should fix the parameters of deformation as in (5.14). So, the 
system of the equations (4.12-4.13) contains the solutions to Bethe ansatz 
equations for both cases. Unfortunetely, we have failed to find an explicit 
correspondence between (4.12-4.13) and (5.15). 



6. Discussion 

In this paper we have obtained the functional relations for eigenvalues of the 
transfer matrix of the si (3) chiral Potts model at q 2 — — 1. In the completely 
homogeneous case we have also developed a direct Bethe ansatz scheme. 

A whole set of functional equations for the usual chiral Potts model has 
been obtained in [26]. However, the si (3) case looks much more difficult. 
In particular, a proper generalization of the Baxter's construction of the Q- 
matrix for the eight-vertex model [5] which works for the usual chiral Potts 
model [6] fails for the sl(3) case. However, the structure of functional equa- 
tions for the arbitrary N should be governed by a proper generalization of 
quadratic "fusion" rules [27, 28] which should have the same functional form 
for the sl(n) case at roots of unity as well. Then it should be possible to 
define a set of Q- matrices for the sl(n) case. Of course, the boundary con- 
ditions for "fusion" rules and analytical properties of solutions at q N = — 1 
will be extremely complicated. However, even a trigonometric limit of the 
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sl(n) chiral Potts model at q N = — 1 is of a great interest, because it corre- 
sponds to the n-layer Zamolodchikov model under the proper modification 
of boundary conditions on a three-dimensional lattice. We hope to address 
these problems in further publications. 
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Figure 1. The square lattice £. 
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Figure 2. The weights iy w (a,/3) and W qp {a,j3) 
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Figure 3. i?-matrix S' a (pi,p 2 ;ri,r 2 ) 




= W™(a,P,>y,6) 



Figure 4. The "star" weight W™Ha, (5, 7, 6) 
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Figure 5. The "star" weight W™'(a, /3,y, 
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Figure 6. The "star-star" relation. 
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